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Abstract 

Studies were performed at INTA in Spain to determine 
water-fdm thickness on a stagnation-point probe in- 
serted in a simulated cloud. The measurements were 
correlated with non-dimensional parameters describing 
the flow and the cloud conditions. Icing scaling tests in 
the NASA Glenn Icing Research Tunnel were then 
conducted using the Ruff scaling method with the scale 
velocity found by matching scale and reference values 
of either the INTA non-dimensional water-fdm thick- 
ness or a Weber number based on that fdm thickness. 
For comparison, tests were also performed using the 
constant drop-size Weber number and the average- 
velocity methods. The reference and scale models were 
both aluminum, 61-cm-span, NACA 0012 airfoil sec- 
tions at 0°AOA. The reference had a 53-cm-chord and 
the scale, 27 cm (Vi size). Both models were mounted 
vertically in the center of the IRT test section. Tests 
covered a freezing fraction range of 0.28 to 1.0. 

Rime ice ( n = 1.0) tests showed the consistency of the 
IRT calibration over a range of velocities. At a freezing 
fraction of 0.76, there was no significant difference in 
the scale ice shapes produced by the different methods. 
For freezing fractions of 0.40, 0.52 and 0.61, somewhat 
better agreement with the reference horn angles was 
typically achieved with the average-velocity and con- 
stant-film thickness methods than when either of the 
two Weber numbers was matched to the reference 
value. At a freezing fraction of 0.28, the four methods 
were judged equal in providing simulations of the refer- 
ence shape. 

Nomenclature 

A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

c Airfoil chord, cm 


*Senior Research Associate, Member 
'Chief, Experimental Aerodynamics 


c p Specific heat of air, cal/g K 
c p.ws Specific heat of water at the surface tempera- 
ture, cal/g K 

d Cylinder diameter or twice the leading-edge 

radius of airfoil, cm 
h Water- film thickness, cm 

h c Convective heat-transfer coefficient, 

cal/sm 2 K 

h G Gas-phase mass-transfer coefficient, g/sm 2 

K Inertia parameter, dimensionless 

K 0 Modified inertia parameter, dimensionless 

LWC Cloud liquid-water content, g/m 3 

M Mach number, dimensionless 

MVD Water droplet median volume diameter, pm 

n Freezing fraction, dimensionless 

p Static pressure, Nt/m 2 

p w Vapor pressure of water in atmosphere, Nt/m 2 
p mv Vapor pressure of water at the icing surface, 
Nt/m 2 

r Recovery factor, dimensionless 

Re Reynolds number of model, dimensionless 

Res Reynolds number of water droplet, dimen- 

sionless 

tf Freezing temperature, °C 

t s Surface temperature, °C 

t st Static temperature, °C 

V Air velocity, m/s 

We Weber number based on droplet size and water 
properties, dimensionless 

We c Weber number based on model size and air 
properties, dimensionless 

We/, Weber number based on water-film thickness 

and water properties, dimensionless 

p 0 Collection efficiency at stagnation line, dimen- 

sionless 

(j) Droplet energy transfer parameter, °C 

A Droplet range, m 

A Stokes Droplet range if Stokes Law applies, m 

Af Latent heat of freezing, cal/g 

A v Latent heat of condensation, cal/g 

p Air viscosity, g/m s 
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9 Air energy transfer parameter, °C 

p Air density, g/m 3 

pi Ice density, g/m 3 

p w Liquid water density, g/m 3 

<J w/a Surface tension of water over air, dyne/cm 

T Accretion time, min 

Introduction 

Since the 1940’s six similarity parameters have been 
identified 1 ' 2 that are important to icing scaling. They 
are the modified inertia parameter, K 0 , defined by Lang- 
muir and Blodgett 3 to insure similarity in droplet trajec- 
tory, the accumulation parameter, A c , to insure 
similarity in quantity of ice accreted, and four parame- 
ters derived from the energy balance: the freezing frac- 
tion, n, the water-energy transfer parameter, 0, the air- 
energy transfer parameter, 9 and the relative heat factor, 
b. The parameters n, 0, 9 and b were defined by Mess- 
inger 4 and Tribus 5 and are not all needed for similarity 
of ice accretion. The AEDC or Ruff scaling method 2 
matches the scale and reference values of the similarity 
parameters K 0 , A c , n, 0and 9. It is applicable to tunnels 
with the capability to simulate altitude. For tunnels 
without pressure control, a modified Ruff method has 
been used in which only the parameters K 0 , A c , n and 
(A arc used. 9 can also be substituted for 0 when this 
leads to more favorable scale test conditions. 

During glaze-ice accretion water does not freeze imme- 
diately in the zone of impingement. The characteristics 
of the resulting layer of liquid surface water are be- 
lieved to influence the shape of the accreted ice to some 
extent and are the subject of a number of ongoing stud- 
ies. Bilanin 0 and Bilanin and Anderson 7 advocated 
adding another similarity parameter, the Weber number, 
to scaling requirements to address surface-water effects. 
Two forms of the Weber number that might be applied 
include that based on the droplet MVD and water prop- 
erties ( We) and one based on model size and air proper- 
ties {We/). A third Weber number was suggested by 
Kind 8 - 910 ; this one was based on the thickness of the 
water film at the leading edge ( We /t ). 

Anderson and Ruff" and Anderson 12 reported the best 
scaling results for glaze ice were achieved when a com- 
promise between constant We and constant Re was used 
with the modified Ruff method. This approach has 
been called the average-velocity method. A possible 
explanation for the success of this method comes from 
the dependence of surface-tension effects on both We 
and Re. For example, studies of droplet-surface im- 
pact 13 for non-icing applications found that splashing 
effects correlated with a K factor depending only on We 
and Reg. Other surface-tension effects such as water- 
film breakup and rivulet formation are also We and Re 
dependent. 


NASA/CR— 2003-211826 2 


Finally, in a recent study at above-freezing tempera- 
tures, Feo and Urdiales 14 measured the water-film 
thickness, h, for heavy-rain conditions. The non- 
dimensional film thickness, h/d, was found to correlate 
with We and Reg The authors suggested that h/d might 
itself serve as an appropriate similarity parameter for 
scaling studies. Although spray characteristics for that 
study were very different from Appendix-C conditions, 
the heavy-rain h/d was evaluated as a similarity pa- 
rameter in studies comparing how methods to select 
scale velocity affected ice shape. 11,12 However, match- 
ing the scale and reference h/d of Feo and Urdiales gave 
scale ice shapes that were generally a poorer match of 
the reference than did other methods of finding scale 
velocity. Later, Feo 15 measured the thickness of the 
water film for MVD and LWC conditions closer to those 
in an icing cloud and established that the new results 
correlated with LWC and Re. With the availability of 
this new correlation for h/d. the constant-water-film 
thickness method needed to be reevaluated. This 
evaluation was a goal of the present study, along with 
the evaluation of the constant-Weber-number based on 
this water-film thickness, We h . 

In this investigation, ice accretion scaling for a model 
size reduction of Vi was evaluated. The Ruff method 
provided the scale test conditions t st , MVD, LWC and T. 
The scale V was found by matching We, We h or h/d to 
the respective reference value or by using the average - V 
method. The study compared how well the scaled ice 
shapes agreed with the reference shapes to evaluate 
each of these methods of determining scale velocity. 

This paper will briefly describe the recent near- 
Appendix C water-film thickness studies at INTA as 
well as the scaling tests performed in the NASA Glenn 
Icing Research Tunnel (IRT). The reference model for 
the scaling tests was an aluminum 61-cm-span, 53-cm- 
chord NACA 0012 airfoil at 0°AOA. The scale model 
was also an aluminum 61-cm-span model but with a 
chord of 27 cm. Both models were mounted vertically 
in the center of the IRT test section. Test conditions 
covered a range of freezing fractions from 0.28 to 1.0. 

Scaling Method 

The basis of the scaling used in this study was the Ruff 2 
method. Ruff demonstrated that the best scaling oc- 
curred when the five similarity parameters, K 0 , A„ n, 0 
and 9, were used. The matching of each of these scale 
parameters to its respective reference value allows the 
solution for scale MVD, time, LWC, temperature and 
pressure. For tunnels without control of test-section 
pressure, either of the parameters 0 ox 9 can be used to 
find scale temperature 12 , with the other parameter ig- 
nored. The choice of scale velocity is left to the user. 



The modified inertia parameter, K n , was defined by 
Langmuir and Blodgett: 3 


K 0 =- + 


A 


8 A a 



( 1 ) 


Equation (6) gives the form of #used by Charpin 1 and 
others. A somewhat more complex form was presented 
by Ruff, 2 and Ruffs expression for 6 was used in the 
calculations in this paper. The values from the two 
expressions are not significantly different, however. 


A match of scale and reference K 0 insures similarity of 
droplet trajectories. In equation (1), AJA Slokes is the 
droplet range parameter, defined as the ratio of actual 
droplet range to that if Stokes drag law for solid- 
spheres applied. It is a function only of the droplet 
Reynolds number, Reg and was tabulated by Langmuir 
and Blodgett. K is the inertia parameter, 

K= P« MVDl V ( 2 ) 

18 dp 

Here d is the diameter for cylindrical models or twice 
the leading-edge radius for airfoils. For the NACA 
0012, a leading-edge radius of 0.0158c was used, where 
c is the model chord. 

The accumulation parameter is: 


Four ways to determine scale velocity were evaluated 
as supplements to the Ruff method. They involved 
either the inclusion of one additional similarity parame- 
ter, or the calculation of scale velocity as the average of 
that required to match scale and reference We and that 
to match Re. The additional parameter was either We, 
We h , or h/d. The We, Wei, and Re are defined as 


We = V2MVD ^ 

(8) 

&wla 


(9) 

^oila 

Re= Vd ? 

(10) 




A. = 


LWCVr 
d p t 


(3) 


The Feo correlation for the parameter h/d will be given 
in the next section. 


A match of scale and reference A c insures that the same 
quantity of ice, relative to airfoil size, is accreted for 
both. This matching assumes that the leading-edge 
collection efficiency, J3 n , for scale and reference situa- 
tions are also equal. However, because J3 0 is a function 
only of K 0 , if K 0 is matched, so too is /?». 

The freezing fraction, from Messinger’s 4 surface energy 
balance, can be written in the form, 

-^K) ,41 


INTA Water-Film-Thickness Measurements 

The water-film-thickness experiments were performed 
by Feo 15 in the INTA 2.8-by-1.7-m wind tunnel shown 
in figure 1. The water spray was generated using a sin- 
gle air-atomizing nozzle located in the open test section 
as seen in figure 1(a). The target probe was positioned 
at the entrance to the bell mouth. Three nozzles were 
used. Spraying Systems models SUJ 12 A, SUJ 12 and 
SUJ 22 B. These operated at different air and water 
pressures and therefore provided some check on the 
validity of the relationship between film-thickness and 


Thus, a match of scale and reference n insures the over- 
all energy balance will agree. The individual terms in 
this expression are (j), the water energy transfer parame- 
ter; 9, the air energy transfer parameter; and b, the rela- 
tive heat factor, introduced by Tribus, et al. 5 These 
parameters are defined as 


<P = tf~t s . 
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2c 


p,ws 


e = \ 


1 1 Lf 
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WW I w 


2c 


A, 


P J 


b = 


LWC V f3 0 c 


(5) 

(6) 
(7) 


- Spray Nozzle 

/ 



(a) Water Spray and Probe. Dimensions are in m. 
Figure 1. INTA Wind Tunnel. 
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cloud conditions. Correlations of the nozzle air and 
water pressures with cloud characteristics at the posi- 
tion of the film-measuring stagnation-point probe were 
obtained by testing with an OAP-200 probe to deter- 
mine MVD and a water-collection probe to measure 
LIVC. The nozzle spray can be seen in figure 1(b). 



(b) Test section with spray. 


Figure 1. Concluded. 

Sixteen tests were performed; the test conditions are 
given in table I. The OAP-200 instrument measured 
droplet sizes no smaller than 25 pm; thus, Appendix-C 
MVD ’ s could not be accurately determined. It was in- 
sensitive to the sprays generated by the nozzles SUJ 12 
and SUJ 12 A, suggesting that MVD’s were below 
25 pm for these. In two cases, the pressures used with 
these nozzles corresponded with the manufacturer’s 
published tests, and these are reported in table I as 
MVD’s of 13 and 15 pm. Test MVD\, then, covered a 
range from well under 25 to 160 pm. 

LWC’s for this study were higher than desired to simu- 
late Appendix-C clouds, but it was hoped that correla- 
tions developed from these data could be extrapolated 
to Appendix-C conditions. The tests were made at am- 
bient temperatures well above freezing, representing a 
freezing fraction of zero. 

A conductance sensor mounted flush in the front face of 
a cylindrical stagnation-point probe was used to deter- 
mine water-fdm thickness. A cross-section sketch of 
the probe 14 is shown in figure 2 with a water-film of 
thickness h indicated conceptually. The probe was lo- 
cated 2.5 m downstream from the water nozzle. A de- 
tailed description of the conductance sensor is given in 
Feo and Urdiales 14 and Feo, Rogles and Urdiales. 16 

Table 1 lists the measured water-film thicknesses re- 
ported in reference 15. The uncertainty in these data is 
estimated to be ±20%. They have been non- 
dimensionalized and plotted in figure 3. Note that 


Table I 

Film-Thickness Test Conditions and Results 


Nozzle 

v, 

m/s 

MVD, 

pm 

LIVC, 

g/m 3 

h, 

pm 

SUJ 12 A 

55 

13 

5.39 

75 

SUJ 12 A 

55 

<25 

4.84 

62 

SUJ 12 A 

30 

<25 

4.53 

79 

SUJ 12 A 

30 

<25 

6.69 

108 

SUJ 12 A 

20 

<25 

6.99 

95 

SUJ 12 A 

20 

<25 

5.17 

89 

| SUJ 12 

55 

15 

5.25 

60 | 

SUJ 12 

55 

<25 

4.22 

65 

SUJ 12 

55 

<25 

2.65 

55 

SUJ 12 

30 

<25 

6.41 

72 

SUJ 12 

20 

<25 

7.79 

106 

SUJ 12 

20 

<25 

8.45 

125 

SUJ 12 

20 

<25 

4.97 

103 

SUJ 22 B 

55 

160 

10.67 

106 

SUJ 22 B 

30 

160 

9.00 

118 

SUJ 22 B 

20 

88 

14.54 

154 

o 


^ d 

L_ 


V 



° 

° 


Figure 2. Stagnation-Point Probe with Water Film. 



{LWC1 pj' 2 Re y \ ! 10' 4 

Figure 3. Film-Thickness Correlation. 
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because of the Re dependency, the slope of the line will 
be affected slightly by test ambient pressure and 
temperature. The best correlation of the experimental 
film thickness data was found to be 


hid 


/ \ 1/2 

' LWC^ 


Pw 


Re- 


in) 


This expression shows no water surface-tension effects, 
indicating that for the spray conditions and the stagna- 
tion probe geometry used, droplet impacts, film break- 
up and rivulet formation are not important. However, 
measurements on the stagnation point on a cylinder or 
airfoil might show a different behavior. 


at 0° AO A. The IRT spray system has rapid-start capa- 
bility; consequently, the model did not require shielding 
during spray initiation. 

In preparing for a test, the temperature and airspeed in 
the test section and the air and water pressures on the 
spray manifolds were set. When these conditions had 
stabilized, the spray nozzle valves were opened to initi- 
ate the spray. The spray was timed for the required 
duration, and then turned off. The fan was brought to a 
full stop and the tunnel entered to record the ice shape. 
A heated ice knife with a cutout in the shape of the 
model was inserted into the leading edge of the ice to 
melt a thin slice down to the model surface. A card- 



(a) NASA Glenn Icing Research Tunnel. Revised Loop with 
New Heat Exchanger Installed in 2000. 


Figure 3 shows that data for all three nozzles 
correlated well with the right-hand side of 
equation (11). The proportionality constant 
that should appear on the right side of equation 
(11) is not known for icing conditions at this 
time. However, it must be some function of 
the freezing fraction, because at a freezing 
fraction of unity, there is no water fdm, while 
at a freezing fraction of zero, the thickness h 
would apply. Furthermore, different constants 
may be required for Appendix-C conditions 
than for the LWC’s of this study. The form of 
equation (11) has been verified analytically by 
Rothmayer 17 for the range of temperatures and 
spray conditions of the INTA tests. 

The INTA tests provided two similarity pa- 
rameters: the non-dimensional water-film thickness, 

h/d, and the Weber number using this film thickness as 
the length scale, We h , defined in equation (9). These 
parameters were evaluated for their use in scaling in the 
IRT scaling tests described next. 

IRT Scaling Test Description 

The scaling tests were performed in the NASA Glenn 
IRT. The IRT is a closed-loop, refrigerated, sea level 
tunnel with a test section of 1.8 by 2.7 m. The tunnel 
loop is shown in figure 4(a). The IRT has 10 spray 
bars. In the summer of 2000, a complete cloud calibra- 
tion 18 for both the mod-1 and standard nozzles was 
completed following the installation of a new heat ex- 
changer and other upgrades. 19 

The models used for these tests were NACA 0012 61- 
cm-span aluminum airfoil sections with chords of 53.3 
and 26.7 cm. They were mounted vertically between 
splitter plates at the center of the IRT test section as 
shown typically in figure 4(b). Horizontal lines at the 
leading edge indicated tunnel center, ±2.5 cm and ±5 
cm from the center as guides for locating ice tracings. 
The mounting arrangement permitted rotation of the 
model for angle of attack changes, but all tests were run 



(b) Model Installed in IRT Test Section. 
Figure 4. Facility and Model Description. 
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board template was placed into this gap and an outline 
of the ice accretion traced. Tracings were taken at the 
vertical center of the tunnel (91 cm from the floor) and 
at 2.5 cm above the center. Ice shape differences be- 
tween the two tracing locations were never significant, 
and only centerline shapes will be reported here. The 
tracings were digitized and the x-y coordinates for each 
ice shape recorded. When plots of ice-shapes were 
compared for different sized models, the coordinates 
were normalized by the model chord. 

Reference test conditions were chosen so that both ref- 
erence and scale conditions would fall within the IRT 
operating envelope. The six cases selected all had the 
same V, MVD, LWC and rbut with temperatures varied 
to provide a range of freezing fractions from 0.28 to 
1.0. Scale conditions for 'A-size scaling were deter- 
mined as discussed above. It turned out that the con- 
stant- IFe and constant- We h methods of finding scale 
velocity gave virtually the same scale test conditions, 
thus, there were only three distinct sets of scale condi- 
tions for each reference case. Nevertheless, tests for 
both the constant- We and constant- We h methods were 
made so that repeatability could be evaluated. Results 
from three test entries — August 2000, November 2000 
and March 2001 — will be presented. Only a sampling 
of typical results will be shown. 

Uncertainty Analysis for IRT 

Tunnel and cloud conditions were recorded every 1 sec 
over the duration of each test for the August and No- 
vember 2000 entries. For the March 2001 tests the re- 
cording interval was 3 sec. Reported conditions are the 
time averages of these records. Estimates of the uncer- 
tainty in the reported average conditions were made by 
considering fluctuations of the values over time, possi- 
ble instrument errors, uncertainties in tunnel calibration 
of MVD and LWC and differences in measurements 
from one location to another in the test section. 

Total temperature was measured as the average of the 
readings of 24 thermocouples located on the D-comer 
turning vanes. Including inherent uncertainty in the 
thermocouples themselves, it is estimated that the total 
uncertainty in temperature is on the order of ±0.5 °C. 

Test-section velocity was determined from the average 
total and average static pressures from 2 pitot-static 
probes, one located on the north wall and one on the 
south wall of the tunnel at the entrance to the test sec- 
tion. The maximum computed difference in velocities 
indicated by the two probes was ±0.25%. Pressures 
also fluctuated modestly with time such that the maxi- 
mum variation in velocity over the spray duration was 
determined to be ±1.3 m/s at an average velocity of 40 
m/s and ±1 m/s for speeds higher than 40 m/s. 


The most significant contribution to LWC uncertainty is 
the calibration process, and scatter in the calibration 
data suggests this uncertainty is about ±10%. 

The largest contribution to drop-size uncertainty comes 
from the instrument used in calibration and was esti- 
mated from the observed randomness recorded during 
calibration. The net uncertainty in MVD is estimated at 
±12%. This uncertainty applies only to the particular 
instruments used (OAP and FSSP probes) in the IRT 
calibration and do not suggest a reference to any abso- 
lute value of MVD, because the absolute value is not 
known. 

These uncertainties in the test parameters produce the 
following approximate uncertainties in the similarity 
parameters for the tests reported here: 17% in K n , 2% 
in /3„, 12% in A c , 15% in n, 1 °C each in <j> and Q, 2% in 
Re, 12% in We, 1% in We h , and 6% in h/d. For the last 
two parameters, uncertainties are in the relative values 
calculated with the correlation of equation (11). 

Scaling Results 

Figures 5 to 11 present the non-dimensional scale and 
reference ice shapes recorded using each method of 
finding scale velocity. Each figure represents a differ- 
ent reference case. Reference ice shapes are shown 
shaded, while the scale shapes are indicated by a solid 
line. The table below each figure gives the test condi- 
tions and similarity parameters for each pair of refer- 
ence and scale tests. The conditions given are the aver- 
age conditions recorded over the duration of each test, 
which can sometimes differ slightly from the planned 
set points. The parameters in the tables were calculated 
from these average conditions. The numbers in the last 
column in the table are not actual water-film thick- 
nesses but rather the relative h/d. This value was calcu- 
lated from the right-hand side of equation (11). 

Scaling Repeatability — Figures 5 and 6 

Figures 5 and 6 each compare two sets of reference and 
scale tests for half-size scaling. Two of these used the 
average - V method of determining scale velocity, and 
two used constant- We or constant- We/,. As noted 
above, these latter methods produce almost identical 
scale test conditions. The freezing fraction was 0.28 for 
the average- V tests (figs. 5(a) and (b)), which were per- 
formed in August and November 2000, and 0.52 for the 
constant -We and constant- We h tests (figs. 6(a) and (b)), 
which were made in November 2000 and March 2001. 

In the August 2000 average - V tests, (fig. 5(a)) the scale 
ice shape closely matched the reference in terms of the 
main ice accretion size and shape as well as the feather 
structure aft of the main shape. When the same refer- 
ence and scale test conditions were repeated in the 
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-0.05 0.00 0.05 0.10 0.15 -0.05 0.00 0.05 0.10 0.15 

X/c x/c 

' 1 Reference Scale 


uncertainties in temperature or 
LWC from run to run can make a 
significant difference in ice shape. 

Large differences in shape are usu- 
ally not seen when the freezing 
fraction is higher than about 0.3. 
Figure 6 shows some typical re- 
peatability results for a nominal 
freezing fraction of 0.5. The refer- 
ence ice shapes shown in fig- 
ures 6(a) and (b) had consistent 
horn angles and ice quantity, as did 


Scaling 
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Figure 5. Ice Shape Repeatability at a Nominal Freezing Fraction of 0.3. 


November 2000 entry (fig. 5(b)), the results showed 
significant differences between the reference and scale 
shapes. In particular, the scale shape appears to have 
been formed at a lower freezing fraction than the refer- 
ence. Figures 5(a) and (b) illustrate how variable ice 
shapes can be for low freezing fraction conditions. At 
freezing fractions as low as 0.28, variations within the 
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the two scale shapes. The results from the November 
2000 (fig. 6(a)) and March 2001 (fig. 6(b)) test entries 
were sufficiently different, however, that the non- 
dimensional reference shape in figure 6(a) was slightly 
larger than the scale, while the scale in figure 6(b) ap- 
peared to be slightly larger than the reference. These 
differences in quantities are consistent with the 
uncertainty in test parameters. 
Thus, when comparing normalized 
scale and reference ice shapes it is 
important to consider test uncer- 
tainties and not expect precise 
matches. 

Rime (Case 460) — Figure 7 

For rime ice, the heat transfer rate 
is high enough that water freezes 
immediately on impact with the 
surface. Thus, for scaling analy- 
ses, heat transfer or surface water 
expressions are unnecessary, and 
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Figure 6. Ice Shape Repeatability at a Nominal Freezing Fraction of 0.5. 
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any method of choosing scale velocity should be 
equally valid. 

Although they serve no purpose in evaluating scaling 
methods, rime tests can be used to validate tunnel cali- 
brations because the different methods of choosing 
scale velocity result in a range of scale velocity and 
L WC. While the scale test conditions are different from 
method to method, they should all produce the same 
rime shape and quantity. 

Figure 7 shows the comparison of scale and reference 
shapes for three methods of determining scale velocity: 
constant We, constant We h and constant h/d. Remember 
that for rime ice h/d = 0, so the values in the h/d column 
in the table have no physical meaning. Nevertheless, 
scale velocities found by matching scale and reference 
relative h/d are valid for use even in rime studies. 

In figures 7(a) and (b) the scale ice shapes simulated the 
reference well with regard to both shape and quantity of 
ice accreted. Even the feather region was well 
matched. The test of the constant-M/ method 
(fig. 7(c)), produced scale accretions with the right 
shapes, but with smaller quantity than the reference. 
Since the uncertainty in LWC has been estimated to be 
±10%, differences in ice quantity from one test to an- 
other should be within about 20%. The differences 
observed in figure 7(c) appear to be of about this 
magnitude. 


Freezing Fraction of 0.76 (Case 461) — Figure 8 

For a freezing fraction of 0.76, only the constant- We 
(fig. 8(a)), constant- Wei, (fig. 8(b)) and constant-M/ 
methods (fig. 8(c)) were tested. All three methods gave 
approximately the same quantity of scale ice and nearly 
the same shape. Feather growth aft of the leading edge 
was well simulated by all methods tested. 

Freezing Fraction of 0.52 (Case 463) — Figure 9 

Shape differences among the different methods are 
more evident when the freezing fraction is reduced suf- 
ficiently to produce glaze horns. Horns were evident in 
tests at a freezing fraction of 0.6. Scaling results at this 
freezing fraction were essentially the same as those at a 
freezing fraction of 0.52, shown in figure 9, with the 
average - V and constant-M/ methods providing scale 
shapes that were good simulations of the reference. 

For n = 0.52, all four methods of determining scale ve- 
locity were tested. Because, as noted above, the con- 
stant- We and constant- We h methods produced virtually 
the same scale test conditions, only the constant- We/, 
results are shown in figure 9(a) for comparison with the 
average - V (fig. 9(b)) and constant-M/ results (fig. 9(c)). 
Neither the constant- We nor the constant- We^ methods 
were able to simulate the glaze horn angle, although 
both produced the correct quantity of ice and repro- 
duced well the feather size and density. 
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Figure 7. Rime Scaling (Case 460). 
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Figure 8. Scaling with n = 0.77 (Case 461). 
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Figure 9. Scaling with n = 0.5 (Case 463). 
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Figure 9(b) shows that again the average - V approach 
provided an excellent match of the horn angle while 
also producing the correct ice quantity and simulating 
the feather region of the reference. Within typical 
shape repeatability, the results for the constant-M/ 
method shown in figure 9(d) were also a very close 
simulation of the reference shape. 


Freezing Fraction of 0.40 (Case 464) 


For this freezing fraction, only the constant -We/,, aver- 
age- V and constant-M/ methods for determining scale 
velocity were evaluated. The best match of the refer- 
ence shape and quantity occurred when scale velocity 
was found by the constant-M/ method (fig. 10(c)). The 
constant- We/, method (10(a)) gave only an approximate 
simulation of the reference shape, while the average - V 
method (fig. 10(b)) produced a poor simulation. Since 
typically the average - V method produces successful 
scaling, the results of figure 10(b) need to be reexam- 
ined with additional tests before conclusions can be 
made. 


Freezing Fraction of 0.28 (Case 465) 


The warmest case tested was case 465, which had a 
freezing fraction of 0.28. Results for three of the four 
methods are given in figure 11. An idea of the variabil- 
ity of ice shapes accreted at freezing fractions this low 
can be seen from the results shown above in 


figures 5(a) and (b). Given the repeatability variations 
demonstrated by those two tests, it can be concluded 
that the constant- Wei, (fig. 11(a)) scale test simulated 
the main ice shape acceptably well. The average- V 
scale results in figure 11(b) do not appear to be a good 
simulation of the reference shape; however, figure 5(a) 
showed results from an earlier test entry for which the 
constant- V scaling simulation at this freezing fraction 
was excellent. The constant -h/d method (fig. 11(c)) 
simulated the main ice shape well, although the scale 
feathers appear to be larger than for the reference test. 

Additional tests at this freezing fraction are needed to 
better assess the constant-M/ method. The preliminary 
assessment is, however, that any of the methods can 
produce acceptable scaling simulation for n = 0.28. If 
scale simulations are needed at such low freezing frac- 
tions, it is recommended that, due to the variability of 
shapes, tests be repeated with more than one test entry. 


Summary 


Measurements of the water-film thickness on a stagna- 
tion probe in cloud conditions simulating those of icing 
led to the non-dimensional thickness expression. 
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Figure 10. Scaling with n = 0.40 (Case 464). 
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Figure 11. Scaling with n = 0.28 (Case 465). 


The expression on the right side of this proportionality 
was used as a relative h/d scaling parameter, and a We- 
ber number based on this relative water-film height, 
Wei,, was also used as a second scaling parameter to 
determine scale velocity to supplement the Ruff scaling 
method. Results from tests with either of these two 
parameters were compared with those from tests with 
constant We and average velocity, which have been 
discussed in previous papers. Scale test conditions 
found using the constant- We and constant- We/, methods 
were the same, within parameter uncertainty. 

Each of the methods typically provided a good simula- 
tion of the size and density of the feather growth behind 
the leading-edge shape for all cases. However, the 
comparison of the main ice shapes near the leading 
edge revealed some differences in how well the various 
methods scaled the ice. The differences depended on 
the freezing fraction. 

For rime ice (n = 1.0) there is no surface water and all 
methods of choosing scale velocity should produce the 
same ice shape. This was generally true, except that for 
the highest scale velocities, the relative quantity of 
scale ice was less than the reference. The deficit ap- 
peared to be within the ±10% uncertainty in the 1RT 
L WC calibration, however. 


At a freezing fraction of 0.76, there was no significant 
difference in the scale ice shapes produced by the dif- 
ferent methods. Scale accretions were somewhat 
smaller than the reference for all methods, but were 
within 10% of the reference size. For freezing fractions 
of 0.40, 0.52 and 0.61, all methods resulted in scale 
accretions closely matching the size of the reference, 
but somewhat better agreement with the reference horn 
angles was typically achieved with the average- V and 
constant-M/ methods than with the constant- We or con- 
stant- Wei,. At a freezing fraction of 0.28, the four 
methods were judged equal in providing simulations of 
the reference shape within the ice-shape repeatability. 
At low freezing fractions, poorer repeatability has been 
observed than for higher freezing fractions. 

Conclusions of this study with regard to the constant- 
We and average- V methods of finding scale velocity 
agreed with those of reference 12, in which those two 
methods were also tested. The conformity of the aver- 
age- V and constant -h/d methods is consistent with the 
presumption that both Re and We are important parame- 
ters in describing those water-film phenomena that af- 
fect ice accretion. 
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